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ABSTRACT 

It is demonstrated that the planar reorientation of a free-free beam in zero gravity space can be 
accomplished by periodically changing the shape of the beam using embedded electromechanical actuators. 
The dynamics which determine the shape of the free-free beam is assumed to be characterized bj the Lu u 
Bernoulli equation, including material damping, with appropriate boundary conditions. The coupling 
between the rigid body motion and the flexible motion is explained using the angular momentum expression 
which includes rotatory inertia and kinematically exact effects. A control scheme is proposed where the 
embedded actuators excite the flexible motion of the beam so that it rotates in the desired sense * w th 
respect to a fixed inertial reference. Relations are derived which relate the average rotation rate 
amplitudes and the frequencies of the periodic actuation signal and the properties of the beam. test 
reorientation maneuvers can be implemented by using feedback control. 

l INTRODUCTION 

Classical models of uniform free-free flexible beams in zero-gravity space result in complete decoupling 
of rigid bod v motion and flexible motion. However, conservation of the angular momentum o 
indicates that the classical models are incomplete in the sense that there is in fact higher order noid - 
ear coupling between the rigid body motion and the flexible motion, if rotatory inertia and kmemati a 1> 
exact modeling effects are included. Assuming that the angular momentum of the beam is always m 
oscillations in the shape of the flexible beam can actually cause a rotation of the beam wi respec 
fixed inertial reference. The rotation of the beam over one period depends only on the shape : of the be., m 
over the period and does not depend on the length of the period; hence thts phenomenon ,s refened 
a geometric phase change. 

T hese observations lead to a scheme for carrying out asymptotic reorientation of a »” 1 * - 

beam in spare rising only electromechanical actuators embedded in the beam. These embeddet 
cha ca Z, nators e g. piezoelectric actuators, do not change the angular momentum of the free-free beam 
buTthev can be "sedfo change the shape of the beam in a periodic way thereby resulting m a rotation ot he 
beam in space. This reorientation scheme, based on the use of embedded actuators, does not requii ■ " 
momentum wheels or gas jets and thus requires a minimal use of fuel to achieve a given beam reonen ,t»i. 

In this paper, the basic modeling issues are addressed. The dynamics which characterize the shape 
of the free-free beam is assumed to be characterized by the Euler-Bernoullt equation, including ma • in 
.lamping, with apppropriate boundary conditions. The coupling between the rigid body motion a, it ■ 
flexible motion is explained using the angulat momentum expression. A control scheme « i proposed u 

the embedded actuators excite the flexible motion of the beam so that it rotates in the desired e. 

Relations are derived which relate the average rotation rate to the amplitudes and the frequencies o 
!L periodic actuation signal and the properties of the beam. These reorientation maneuvers - 
implemented by using feedback control. Important features of the approach are indicated. 


2. A PLANAR FREE-FREE BEAM MODEL 


Consider a uniform free-free beam of undeformed length 2 L in space with zero angular momentum and 
zero linear momentum. Referring to Fig. I the motion of the beam is constrained to a plane defined by 
vectors (61,63) where (et^.ta) is an orthonormal basis for an inertial frame whose origin is at the center 
of mass of the beam. Let (i*j,k) be a rotating frame with its origin fixed at the origin of the inertial 
frame such that the vectors (LA:) lie in the plane ( e i . ^3 ) and j = 62- The straight line passing through 
the origin in the direction of vector A* is called the reference line. Let the beam initially be at rest in a 
straight line configuration aligned with the reference line. Then, the location of each point on the line of 
mass centroids of the beam can be described in terms of the parameter s € [—L.L\. This parameter $ 
can be viewed as a label for each of the rrossections. We assume that as the beam deforms the shape and 
the area of the crossections remain invariant. Following other researchers 1,2,3 we introduce three functions 
y(s, t) : [-L, L] x ft — ft and v( s. t ) : [-L, L] x ft — T x such that (u(sj) + .s, y(sj)) define t he 
coordinates of the line of centroids in the deformed configuration with respect to the moving frame (L A) at 
time t. The angle w(sj) between the normal to the crossection at s and £3 specifies the orientation of the 
crossection. The normal to the crossection at s is denoted by t 3 . We define the material basis (^,^2^3) to 
be orthonormal so that t\ lies in the plane (61,63). The crossection itself can be associated with the set of 
points (£1, £2) in a compact set A C ft 2 such that fi t\ + £2^2 + (u(.s, t) + s)k + (y(s y t))i gives the location 
of any point on the beam as £1 and £2 vary through A and s varies from —L to L. 


Since the origin of the inertial frame is fixed at the center of mass of the beam we obtain 


J y(s,t)ds = 0, 
J u(s,t)ds = 0. 


(U 

( 2 ) 


Let, p denote the constant mass density per unit volume of the beam. We assume that the beam has a 
symmetric crossection so that the first moment of inertia of the crossection about the line of centroids is 


J Ptidhdh = 0 . 

The second moment of inertia of the crossection about the line of centroids is 

h = J rfMxdh. 


The mass per unit length of the crossection is given by 


m 0 - 


j pd£\d£ 2 . 


O) 


(')) 


We define the angle 6 {t) between e3 and k so that y(s,t) measured from the reference line satisfies tin 


following orthogonality condition 




sy(s, t)ds = 0. 


! f) ) 


The existence of the angle 9 (t) follows from the geometry indicated in Fig. 1. This definition provides a 
separation between the motion which determines the shape of the beam, given by y(s,t ), - L < s < L. and 
the rotation of the beam as a whole, given by 6 (t). 
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3. EQUATIONS OF MOTION 


We first develop a kinematically exact expression for the angular momentum of the free free beam. Let 
6i.0. t) he vector from the origin of the inertial frame to a point on the beam at time 

t; then 

j = (s Sin 0 + y cos 9 + cos * + u sin 9)e x + (6)*2 + (s cos 9 - 6 sin 0 - y sin 9 + u cos 9)e 3 (7) 

where 9 = 9(t), y = y(sJ) and vf = The angular momentum about the origin of the inertial frame 

at time t is zero so that 

J L J^ x d -g d(l d( 7 is = o. m 

Substituting equation (7) into equation (S) and using equations (4) and (5) we can express » in terms of 
i). u and a as - a , . 

f^rA m o^ + [it + m °^ u ~ 


9 = 


jt L {-moS 2 - m 0 y 2 - h}ds 


(9) 


where a = v - 9 is the angle between the normal t 3 to the crossection at s and the reference line. 

Assume that the beam is unshearable and inextensible and that the deformations are small. This 
implies, using equation (2), that 

u(s, t) = 0. 1 

and that / 1 i \ 

a^y s . 1 

We use the Euler-Bernoulli beam model to characterize the shape of the beam. Thus y{s,t) satisfies 
Euler- Bernoulli equation of the form 


moytt + 7 Vtssss + EIy S9 9s - ( 5 

j= i 


with the boundary conditions 


yss(~L) — VssiL) — 0 , 

J Jsss{ — L) = ysssiL) — 0 


( 1 * 2 ) 

< i;r. 
i\\) 


where I = I 2 /p, E is Young’s elasticity modulus, S' is the distributional derivative of the delta function and 
where for simplicity we assume Kelvin- Voigt damping with a positive damping coefficient 7 - In ^dd, ion 
y(s.t) must satisfy conditions (1) and (6). Internal bending torques t b \t), j = 1,. • • ^ ™ ^ 

m point actuators located at , = on the beam where s, € [-L,L). These ^dedde^ 
actuators change the shape of the beam but at the same time preserve the angular momentum. A1 
such actuators are capable of inducing relatively small displacements one can excite the beam period all 
at a frequency near one of the lower resonant frequencies of the beam to obtain relatively g 1 

shape change. 

Using expressions (6), (10) and (11) in equation (9) we obtain 


9 = 


~ f-L hytads 
r + fti m 0 y 2 ds 


i 1 > ) 


where r = §m 0 I 3 + 2I 2 L. This expression demonstrates the nonlinear coupling between the beam's shape 
and its rigid body motion. Expression ( 15) is non-integrable in the sense that if y(s.t) is a periodic function 
of time, the integral of 9 over one period is, in general, non-zero. 


Remark 3.1 If in the above derivation we had not used the kinematically exact expression for the angular 
momentum but had used the linearized strain assumptions we would have obtained the expression 


1 f L 

9 = -- hlhtds . 
T J-L 


( 16 ) 


As can be seen expression ( 16) leads to the incorrect conclusion that a periodic change in the shape of the 
beam does not result in rotation of the beam. Note that inclusion of rotatory inertia effects and the use 
of the kinematically exact expression for the angular momentum is necessary in order to demonstrate that 
the beam can rotate in space due to periodic shape change. 


We expand the solution y{sJ) to equation (12) in the series 


oo 

y{s,t) = ^2w t (s)qi(t) U") 

: = 1 

where Wi{s)* i — 1,2, . . . are the orthonormal elastic mode shapes of the Euler-Bernoulli model. The elastic 
mode shapes are given by 


Wi(s) 


COS( fiiS) — 

sin (&{$) + 


sin(/3, L) 
sinh(/3,L) 
cos{j3 t L) 
cosh(/3, L) 


cosh {j3is) 
sinh(/?t^) 


if i = 1,3,5... 
if i = 2,4,6... 


where J, are the positive roots of the equation 

cos('2t3L ) cos\\{'23L) — 1 


ordered according to their magnitude. 

Expansion (17) provides the modal decomposition 

m 

Qi T “b Q = 'y " j Uj (i ), t — 1 , 2 . . . \ 1 ^ ) 

where c t = and b tJ = Equation (12), or equivalently equation (18), determine 

the shape of the beam and is called the shape space equation. Substituting equation (17) into equation i 15) 
we obtain 


A zhTM Ml i ioi 

r + Eg ,?. 2 

where = w t (L) - Wi(-L). We note that (19) is non-integrable for any truncation of the infinite series 
in ( 17). 
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4 ASYMPTOTIC reorientation maneuvers 


The goal is to accomplish asvmptotic maneuvers, i.e. starting with Wo) - yM) - yd».*o) - 0 
we want to rotate the beam to that 0(f) - 0 4 , !/(*,<) - 0 and *(».() - 0 as ( - oc for some ties, red 

angle tf.v. 


Consider the periodic excitation of the beam at a single frequency u as 

i-j(t) = v° + v'j cos(u;t), j = 1,2. ...,m 


( 20 ) 


Since the shape space dynamics of the free-free beam is asymptotically stable, the steady-state motion of 
the beam is given by 

q t (t) = /, + a, cos (wt + <M 


where 


', = ^7 EM- 

• J=1 


( 21 ) 

( 22 ) 


and 


a, = 

0W, - w 2 ) 2 + c, 2 u;, 2 j=i 


. , / <W 

— CLTCIQ 1 9 2 

1 ' L? • LiJ 


(23) 


(241 


, 

The excitation function (20) should be sufficiently small so that the Euler-Bernoulli model for the shape 
space dynamics remains valid 6 . 


If qf is small comparing with r we can approximate 


and thus 


r + ESi 9? ’’ 


1 00 . 1 
$ ~ " '^1 d* 2 

z = l 






Li=l 


£C 


i=l 


Integrating over one period and using equation (21) we obtain 

dt 

4>x) 


1 ^ 

1 ^ 

1 

0 

II 

u- 

00 

OO 

X>9. 

EC 

^ JO 7 

.t=l 

J = X 


r— 1 00 

= / -5 y'-a,J,u;sin(wt + 

r L.=i 

O- 00 00 

= — jT aiJ,lj a jSin(<l>j - <t>i)- 


-3 

E<'> + a ; COS(lJ* + 0 ; )) 
j=l 




Expression (25) implies that, in general, the change in angle 0 in steady-state over one period is non-zero 
thereby proving that a periodic change in shape of the beam results in a rotation of the beam; the steady ^ 
state difference 0(^j) ~ 0(0) is referred to as the geometric phase. There are cases, however, w len ic 



geometric phase turns out to be zero. 


Proposition 4.1 Assume that the steady-state motion of the beam is described by equation (21). Then* 
8{ — ) - 0(0) = 0 if any of the following conditions hold: 

1. a t — 0 for all i 

2. /, = 0 for all i 

3. o t = ci>^ for all i , j 

The second statement of the proposition is the most important. It implies that for a non-zero geomet- 
ric phase the beam should necessarily vibrate about a non-straight line reference configuration. It follows 
from expression (19) that following the motion q t (t) = — a,- cos(cjf + <&) yields a steady-state geometric 

phase change negative to that of (21). Therefore, in order to rotate the beam in the opposite direction it 
is sufficient to reverse the signs of v‘f and 

Remark 4.1 Expression (25) can be used in order to predict the sign and the value of the geometric 
phase. Consider a beam which has a square crossection with side size R. Assume that two actuators at 
Sj = -rL and = rL where 0 < r < 1 produce torques according to equation (20). Using two first modes 
in the series (25) yields 

0 ( — ) 9 ( 0 ) - 81.5661((t>r) 2 - (v?) 2 )K - t;g)S(r)sin(<k - fr) j2|j) 

w ' p 2 L' i ER 7 {\{ £) 2 + \ - u ; 2 ) 2 + c\v 2 yj(ul - W 2 ) 2 + C^u > 2 


where 


<f>(r) = (sin(2.36502r ) 4- 0.1329 sinh(2.36502r)) 2 (cos(3.9266r) - 0.0279 cosh(3.9266r )). 

We are now in a position to formulate a specific control strategy to accomplish the asymptotic ma- 
neuver. Starting at rest with 8{t 0 ) — 0o application of control law (20) results in a nonzero rotation nvor 
a period. By repetition of cycles of motion (21) as many times as necessary the beam can be caused to 
rotate closer and closer to 84. As 0(f) approach 8^ we can reduce the amplitude of the oscillations to zero 
in a way so that 8(t) — 8^ as t — ► oo. 


The proposed control law is of the form 


Vj(t) = 5fc 


+ v'j cos(ut) 


= 1, 


m. 


where 2{k < / - f 0 < k - 1,2,...; that is, the control excitation is an amplitude modulated 

function, where vj are constants and Sk denotes the scalar amplitude modulation sequence that dmin**> 
the control excitation on the k- th cycle. Each cycle is exactly p periods. 


The constants u;, v*j , v'j can be chosen nearly arbitrary, although one approach is to choose r u . r* to 
maximize the geometric phase expression 

CO 'X) 

X ] ^2 a, J.ljOj sm(4>j - <?,) 
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Since |- fc | - 0 then q t (t) - 0 and q, - 0 as t - oo. By continuity 9(t) - 6 con for some constant 0"" 
as t — o c. We want to show that 0 con = 0,(. 


By contradiction, assume that 6<° n > 9 d . Let h > 0 be sufficiently small so that 9™ n ~h > Choose 
so that - c x 

>g>o. 

A0* / 


Then, there exists an integer such that for any k > N 3 it follows that |f fc | < 6 and - 0 eon | < h- 
Note that for any k > iV 3 + 1 and I > Yj + 1 


'•fc = 


ft* ~ ^ e , 

Ad' 


< 0 


and 


r’t 


> 


h~9 d 


A0* 


)‘ 


> > k/l 


Thus, we conclude from (A2) that for any k,l > tV 3 + 1 it follows that £fc £l # 0- Hence, we obtain 
a contradiction to the convergence of the sequence e k to zero as fc - oo. Similar arguments lead o a 

contradiction in case 9 LDn < 9,(. ^ 


Finally, it follows from equations (28) and (23) that 


lim 6{ t ) = Od* 

t — 50 t-+rx> 




= 0 < $ < L 


The controller which we have constructed has two functions. Its main function is to excite the om illa- 
tions of the beam in such a way so that the beam rotates in the desired sense. Subsequently, the cont.olhu 
serves to suppress the vibrations previously excited so that the free-free beam comes to rest with a desi.ed 
orientation. Note that control law (27) is a non-smooth feedback control law. 

A NUMERICAL EXAMPLE 

Space structures can often be modeled as light and flexible beams. Consider a beam with half-let.-o h 
L = If ml, density per unit volume p = 1400[/rg/m 3 ] and square crossection with the side size ft - 0.1,^.. 
Young’s modulus of the beam is E = 3.0 X 10 6 [A/m 2 ] and the damping coefficient of Kelvin- Voigt dampm* 
is -v = 0 2. Two actuators are installed near both ends of the beam at r = 0.9. The maximal torque 'Mill 
of the actuators can produce is equal 100[rVm]. The excitation frequency w = 13[//z] is selected to he 

_ - . _ „ r ^ 1 . r i— fi u onH J — [ , J ,1 1 C 


OI Uie aCLUcUUl& CdU ^ W L * J * _q * 

between the first 10.6[flz] and the second 29 [Hz] resonant frequencies of the beam; v } and 
chosen using expression (26) to maximize the geometric phase change over one period. For this examp - 
we choose p = 5 and 7l = 7 2 = 0.9. We want to rotate the beam from 6 0 = 0A[rad) at I t - 0 [>" *> 

6 d = Ofradl. The dependence of the angle 9(t ) [rad] on time <[sec] is shown for a part of the maneuvm m 
Fig. 2. In this case the geometric phase change over one period in steady-state predicted by expressio 

o ° . t 1 • 14 .^ O A 1 1 1 w 1 \T*n H\ T IP npnp 


tie, L. ill LIUS case Ulc gwiucuiL r - 4 r n mi J 

is e„„al to -2.7465 x 10- [ruti] whereas its actual value is equal to -3.0411 X 10- [rad). The del,™, 
of the control parameter s on time is shown in Figure 3. 


Ml \ J * ) 
(Mire 


where a,, /,,</>,, t = 1, . . are related to vf, v*,j = 1 m according to expressions (22)— (24). and r°. r* 

are constrained by 

m m 

j=i j=i 

In terms of Vj % = 1, m this is a constrained mathematical programming problem which is linear 

in Vj (for fixed vf ) and quadratic in (for fixed v®). We will subsequently denote the maximum value of 
this constrained optimization problem as A9* . 


The modulation sequence £;.+ 1 


is defined in terms of an "average” of 9{t ), over the A;-th cycle, that is 


9' l k ve = -(ma xO(t) + min 9{t)) 


(2-S ) 


where the maximum and minimum are over 2 ^ k 1 - ?f " < t — tn < 'QUL m 

jj — — jj 

variables 9q vc = #o and £o = sign ( 9 ) • We express Sk in terms of 9 k v J x 
(Al) Compute 


r k = 




We also introduce two auxilary 
and c 1 as indicated below: 


(A2) In case |r*| > if and have the same signs then £*; = kfc-i |sign( ) ; if r* and have 

opposite signs then £* = 7ikfc-i|sign(rfe), where 0 < 71 < 1. 

(A3) If 0 < \ r k\ < kfc-il then Sk — l 2 r k< where 0 < 72 < T 

(A4) If r k = 0 then Sk = £^_i . 


Proposition 4.2 If the proposed control law is of the form (27) where Sk 
steps (A1)-(A4), then 


lim 9l ve — 9d, lim Sk — 0. 

k — oo k—>o o 


is selected according to 


Sketch of the Proof. The sequence \sk\ is non-increasing and bounded on [0. 1]. Therefore, there exists 
6 G [0. 1] such that b = inf ^ \sk I- We want to show that 6 = 0. 


By contradiction, assume 6^0. Then, for £ = ^2^1^) ^ we can an i ntG S er ^1 suc ' 1 r ^ iar * or 
all k > .V, Ifffcl - b < £. From (A2) and (A3) we conclude that only two cases are possible: Sk = 6 for all 
k > Ai or Sk = “6 for all k > AV 

Assume that the former case is true. Since the transient decays to zero and using continuity of 9 with 
respect to q t and q t we assert that for £1 = ^6 3 A0* there exists an integer N 2 such that for any k > V_> 

b 3 ±e’ - h < 9 a k i\ - or < b 3 ±e m + fi, 


'--•SrV 


where A0* > 0. Note that > 6 3 > 0. Choosing an integer / so that / > 2 6 ^ -f- + 1 we 


com 


huh* 


that 


0 > 


)/ + *«,- ^ ^ e d - e^ +i ^ u3 

A6>* > A 6' 


> b 3 > 0. 


Therefore, the former case can never occur. Similarly, we can verify that the latter case also leads r< 
contradiction. Hence, 6 = 0. 
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6 CONCLUSION 


In this paper the angular momentum expression for a planar free-free beam in space is derived It is 
shown how the general motion of the beam can be separated into rigid and elastic motions. The change 
of shape of the beam is described by the Euler-Bernoulli equation with free-free boundary conditions 
Angular momentum conservation leads to the nonlinear dependence of the rigid motion on the shape of 
the beam As shown this dependence is non-integrable in the sense that a periodic change in shape of the 
beam results in a non-zero rotation of the beam over one period. Approximate relationships expressing the 
average rate of rotation of the beam in terms of the amplitudes and phases of periodic excitation of the 
beam by internal actuators are derived. Finally, a control strategy for a planar asymptotic reorientation 
maneuver is developed. 

A general treatment of the interplay between deformations and rotations of deformable bodies is given 
bv Shapere and Wilczek. 7 Reyhanoglu and McClamroch 8 have developed a framework for reorientation 
of multibody systems in space. In this paper, we have used the framework developed by Shapere and 
Wilczek for the specific problem of reorientation of a free-free beam in space; our results represent, in a 
certain sense, the limiting case of the multibody results obtained by Reyhanoglu and McClamroch w en 
the number of bodies increases without limit. 

Although our study in this paper has been concerned with the ideal case of reorientation of a free-free 
beam in space, we note that the same ideas are applicable to reorientation of a wide class of deformable space 
structures, using only actuators embedded into the structure. In this sense, smart structures techno ogy 
can be used to accomplish a variety of efficient reorientation maneuvers for space structures. 
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